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rn ■ ABSTRACT 

o 
o 

I We have conducted a deep multi-color imaging survey of 0.2 degrees^ centered on 

' the Hubble Deep Field North (HDF-N). We shall refer to this region as the Hawaii- 

Q ■ HDF-N. Deep data were collected in U, B, V, R, I, and z' bands over the central 0.2 

' degrees^ and in HK' over a smaller region covering the Chandra Deep Field North 

. (CDF-N). The data were reduced to have accurate relative photometry and astrometry 
across the entire field to facilitate photometric redshifts and spectroscopic followup. 

^ ! We have compiled a catalog of 48,858 objects in the central 0.2 degrees^ detected at 5a 

I significance in a 3" aperture in either R or z' band. Number counts and color-magnitude 

I diagrams are presented and shown to be consistent with previous observations. Using 
color selection we have measured the density of objects at 3 < z < 7. Our multi-color 

■ data indicates that samples selected at z > 5.5 using the Lyman break technique suffer 

■ from more contamination by low redshift objects than suggested by previous studies. 
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Subject headings: galaxies:observations, galaxies:evolution, cosmology:observations 

1. Introduction 

Deep surveys provide numerous constraints on the structure and evolution of the universe. It 
has been known that galaxy number counts in a specific bandpass can provide a useful constraint 
on galaxy evolution (Tinsley 1972; Brown & Tinsley 1974; Tinsley 1980). Multi-color surveys 
further constrain the galaxy formation history and allow for the selection of galaxies in specific 
redshift ranges. This work was pioneered in the Hawaii surveys (Cowie 1988) and subsequently 
used extensively by Steidel et. al (Steidel et al. 1996, 1999) for mapping z ~ 3 and z ~ 4 galaxies. 
By using a color selection these groups were able to select and obtain spectra for over 1000 z ~ 3 
galaxies with a 90% success rate. This gave us a wealth of information about this redshift range. 
With the ability of modern detectors to take deep images in the near IR we can now select galaxies 
up to z > 6.5 (Hu et al. 2002). The Hubble Deep Fields (HDFs) (Wilhams et al. 1996) showed the 
value of photometric redshifts applied to deep multi-color surveys. This technique allows estimates 
of redshifts for objects much too faint for optical spectroscopy. It allows for economical measurement 
of the redshifts of millions of galaxies. 

Despite the depth and accuracy of the HDF data, it only covers a small co-moving volume. 
Postman et al. (1998) have looked at the clustering of galaxies using the two point correlation 
function of galaxies projected on the sky. This work covered 16 sq. degrees in / band. It has shown 
that galaxies are highly clustered on scales of 15-20 h^^ Mpc where h = Ho/100 km s~^Mpc~^. 
This has provided a challenge in obtaining an unbiased sample of galaxies, since several square 
degrees must be surveyed. The present paper is the first of a series which will aim to obtain the 
required data by imaging several 0.2 sq. degree fields. 

Previous studies of galaxy formation and evolution have solely used optical data which intro- 
duced biases. However new generations of radio. X-ray, and mid-IR telescopes have opened up 
these wavelengths to deep wide-field surveys. X-rays have constrained the Active Galactic Nu- 
clei (AGN) history (Barger et al. 2002), while radio and mid-IR fiuxes have proven to be a good 
probe of star formation. Space based optical imaging also allowed for the morphological study of 
galaxies out to high redshift. As such we have chosen to target regions of the sky that have been 
observed in multiple wave bands. We began with the HDF-N because of its deep and complete 
multi-wavelength data provided by the Great Observatories Origins Deep Survey (GOODS)^^ and 
other surveys. We have covered 0.4 square degrees centered on the HDF-N. The central 0.2 square 
degrees of this data are of consistent quality in all optical bands, this area contains 48,858 objects 
detected at 5a in R or z' band. 

This paper is one in a series focusing on the Hawaii HDF-N. The properties of X-ray selected 
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objects are discussed in Barger et al. (2003, 2002). Photometric redshifts, spectroscopic redshifts, 
and multi- wavelength analysis will be discussed in Capak et al. (2003). The main focus of the 
present paper is the optical/IR data reduction and catalog. We also provide number counts and 
color-color plots to allow for comparison with current and future surveys. We have also selected 
high redshift candidates using color selections and comment on the effectiveness of this technique 
at 2; > 5 as well as the implications for the star formation history. 

2. Data Reduction 

To optimize our data collection, observations were conducted on a variety of instruments. The 
U band data were collected using the Kitt Peak National Observatory 4m (KPN0-4ni) telescope 
with the MOSAIC prime focus camera. This camera has a reasonable U band response and 36' x 36' 
field of view (Jacoby et al. 1998; MuUcr et al. 1998; Wolfe et al. 1998). The B, V, R, I, and z' band 
data were collected using the Subaru 8.2m telescope and Suprime-Cam instrument (Miyazaki et al. 
2002) which is optimized for red optical response and has a 34' x 27' field of view. Finally our HK' 
data were collected using the QUIRC camera on the University of Hawaii 2.2m telescope (Hodapp 
et al. 1996) with a 3.6' x 3.6' field of view. The HK' filter covers both the H and K' bands in a 
single filter which allows for greater depth at the expense of some color information. The details 
of these observations arc found in tables 1 and 2. 

For data collected with Suprime-Cam a five point X-shaped dither pattern was used with 
one arc minute steps. The camera was rotated by 90 degrees between dither patterns to remove 
bleeding from bright stars and provide better photometric calibration. The KPN0-4m data were 
collected using a nine point grid as a dither pattern with one arc minute steps between pointings. 
For both the Suprime-Cam and KPN0-4m observations the telescope was offset by 10" in a random 
direction between dither patterns. The HK' data were collected using a 13 point diamond shaped 
dither pattern with 10" steps between exposures. The filter profiles multiplied by the detector 
response are plotted in Figure 1 along with the HDF-N filters used in Fernandez- Soto et al. (1999). 
Quantifiable measures of image quality are given in Section 2.4. 

2.1. Flat Fielding 

All of these data were reduced using Nick Kaiser's Imcat tools^"^. The U, B, V, R, /, and z' 
data were first overscan corrected and bias subtracted. Median sky fiats were then used to provide 
a first pass flat fielding. Objects were then masked out and a second fiat field was generated 
for each dither pattern. The second pass flat fielding was necessary to correct for changes in the 
instruments due to mechanical flexure and rotation of the camera and optics. The sky was then 
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subtracted using a second order polynomial surface fit to the sky for each chip of each image. After 

the first pass sky subtraction the Suprime-Cam images had time variable scattered light structure 
left in the images. This problem occurred at the edges of the field where the camera was vignetted 
and was worst when the moon was up. To correct for this structure the objects were masked out 
and a surface was tessellated to the sky background and subtracted off. 

The z' band data taken before April of 2001 were corrected for fringing. To achieve this a 
surface was tessellated on a 32-pixel grid over the first pass flat field. This surface was used to 
provide first pass fiat fielding instead of the median fiat. A median fringe frame was then generated 
for each night of data, scaled to the background in each image, and subtracted from the fiat-fielded 
images. The images were then second pass flattened and sky subtracted as with the other bands. 
After April of 2001 new chips were installed that had minimal fringing so fringe subtraction was 
not performed. 

The U band data were corrected for internal pupil reflections which occur on the MOSAIC 
camera on the KPN0-4m telescope. A pupil image was constructed by dividing the U band flat by 
a V band flat and masking out the areas that did not contain a pupil reflection. The pupil image 
was subtracted from the IJ band flat before flat fielding. The pupil image was then scaled to the 
background in each image and subtracted off. A second pass fiat fielding was then performed for 
each dither pattern as done for the other bands. 

The B.K' data were collected in 13 point dither patterns so the data could be fiattened and 
sky subtracted separately for each dither pattern using median sky flats. This allowed us to remove 
the time variability in the bias, flat field, and sky which occurs in IR arrays. The images in each 
dither pattern were then combined using a weighted mean to provide an image deep enough to map 
onto the optical astrometry grid. Cosmic rays were removed during the combination process by 
looking for pixels which lay more than 5o" from the weighted mean value. Sigma was calculated from 
the background noise in each image. Since the image depth varied across the combined, dithered 
images, an inverse variance map was created for each dither sequence. These inverse variance maps 
were later used to combine the separate dither sequences. 

2.2. Astrometry 

An initial solution for the optical distortion was calculated by comparing standard star frames 
to the USNO-A2.0 (Monet et al. 1998). From this point onward we treated each chip in each 
exposure as a separate image with a separate distortion. This took into account changes in the 
optical distortion with time, as well as any mechanical motion on the focal plane. 

Using the initial solution we registered images to one another to cross identify point sources 
from image to image. We also identified any USNO-A2.0 (Monet et al. 1998) stars which were 
not saturated in our images. Since most USNO-A2.0 objects were saturated in the Suprime-Cam 
images, the MOSAIC U band data were used for the astrometric solution. The IJ band data was 



-5- 



also chosen because it covers a larger area than the other bands. 

At magnitudes fainter than 18 the USNO-A2.0 stars are not properly corrected for magnitude 
dependent systematics. As well errors for individual stars are not quoted in the USNO-A2.0 and are 

significantly worse than the mean at the faintest level. As such all USNO-A2.0 stars were assumed 
to have a Root Mean Square (RMS) uncertainty of 2" in relative astrometry which is larger than 
any expected systematics. 

In the center of the field USNO-A2.0 stars were only used to constrain the scale factor and 
absolute astrometry of the fit. The relative astrometry was calculated by minimizing the positional 
scatter of stars which appeared in multiple images. This was done by fitting a third order two- 
dimensional polynomial to each image. The coefficients of the polynomials were calculated by 
minimizing the chi-squared defined by equation 1 with respect to the coefficients defined in 
the polynomial function. In equation 1 iVj is the number of images, and Ng is the number of stars 
used for the astrometry. In this equation the USNO-A2.0 is treated as an image, however f{x,y) 
is replaced with the known star positions. This effectively constrains the scale factor and rotation 
of the image. Stars which were incorrectly cross-identified between images were removed in an 
iterative fashion until the astrometric solutions converged. 



The resulting relative astrometry should be good to the centroiding error of ~ 0.1 pixel or 
« 0.03" RMS across the center of the field where there are many measurements for each source. 
However this method breaks down at the edge of the field where there are fewer images. At the 
very edge of the field the astrometry is only good to 0.5" since only the USNO-A2.0 constrains the 
fit hence systematics can not be removed. 

Any chromatic aberrations in the astrometry should be minimal since both the KPN0-4m and 
Suprime-Cam have an atmospheric dispersion corrector. A large number of objects, with a wide 
range of colors were used for the astrometry, so any residual effects should be averaged out. A more 
extensive discussion of this issue can be found in Kaiser (2000). 

Once an astrometric grid of objects was established, the B, V, and, R data were warped onto it 
using a two-dimensional, third order polynomials for each image (still treating each chip separately). 
To improve the grid in the redder bands point sources from the astrometricaly calibrated R band 
image were added to the U band astrometric grid. The /, z', and HK' data were warped onto this 
improved grid in the same way as the bluer data. To provide accurate absolute astrometry the final 
images were registered to the radio catalog of Richards (2000) . 

We have included comparisons to the radio catalog and USNO-Bl.O (Monet et al. 2003) in 
Figures 2, and 3. The USNO-Bl.O positions are systematically offset from the radio positions by 
0.5" north and 0.2" east (see Figure 3), but this offset was removed before making other comparisons. 
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The RMS astrometric scatter between the our positions and the radio positions over the central 0.2 
square degree field is 0.22" . It is 0.32" between our positions and the USNO-Bl.O. In both cases 
the scatter is dominated by the astrometric errors in the reference catalogs which are 0.16" in the 
radio catalog and 0.34" in the USNO-Bl.O and shows no systematics across the central 0.2 square 
degree field. As expected we observe systematics with respect to the USNO-Bl.O catalog at the 
edges of the field, however these offsets are less than 0.5" in amplitude. Furthermore systematics 
at this level are known to exist in the USNO-Bl.O (Monet et al. 2003) so we chose not to remove 
them. 



2.3. Image Warping And Combination 

For the U, B, V, R, I, and z' data chip-to-chip photometric scaling factors were calculated 
for each dither pattern. This was necessary because each dither pattern was flattened separately 
which can change the relative normalization of the chips. When collecting our data we used large 
steps in our dither patterns and rotated the camera when possible. This meant we could always 
find a reference chip which covered the same area of sky as two or more chips in another exposure. 
The scaling between two chips in one exposure could then be calculated by comparing their scaling 
factors to a chip in another exposure. This was independent of exposure to exposure scaling so 
long as the photometry did not vary across the reference chip. To calculate the scale factors we 
used 6" apertures on objects detected at greater than 50a with half light radii less than 1.5" . 
The large apertures were used to prevent any variations in seeing from affecting the scaling. The 
large number of exposures meant we made many measurements of the chip-to-chip scaling between 
adjacent chips. These were combined in a weighted mean. Once the chip-to-chip photometric 
scaling was removed, exposure to exposure scale factors were calculated in a similar fashion by 
comparing overlapping areas between cxposTircs. 

The geometric distortions in the MOSAIC and Suprime Cam optics cause each pixel to have a 
different effective area on the sky. By flat fielding these images we introduce a geometry dependent 
photometric offset. When calculating the photometric scaling we corrected for this effect. 

The HK' data consists of 210 dither patterns of 13 images each. These data were collected 
under a combination of photometric and non-photometric conditions. We collected these data 
in such a way that each dither overlapped several neighboring dither patterns. To correct the 
photometry we calculated scale factors for each dither pattern such that the photometric scatter 
was minimized across the field in the overlapping regions. The absolute zero point was allowed to 
float during this procedure and later tied to the Fernandez-Soto et al. (1999) catalog (see §2.5). 

Saturated pixels were clipped in the U band images before warping or combination. However 
saturated pixels in the Suprime-Cam behaved in a very non-linear fashion, often dropping or ringing 
in value as they became more saturated. As such we could not simply clip the saturated pixels. 
We removed the saturation spikes along with satellite trails by searching for strings of elongated 
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objects which formed straight hnes in the images. A strip 21 pixels wide was cUpped from the 
images along detected lines. 

After applying photometric corrections each image was warped onto a stereographic projection. 
We corrected the geometry dependent photometric scaling by choosing a mapping which preserved 
surface brightness. The re-sampling was done using a nearest neighbor algorithm. A weighted mean 
was then calculated for each pixel, pixels which lay more than 5a from the mean were rejected to 
remove cosmic rays. For the U, B, V, R, I, and z' data the sigma used in the rejection and weighting 
was measured from the RMS background noise in the un-warped image. For the HK' images a 
sigma was calculated for each pixel during the dither pattern combination and carried through 
to the final image combination. To avoid clipping pixels which varied due to seeing variations a 
window was defined around the median of the pixel values. If a pixel value were between half to 
twice the median value of all pixels it was accepted even if it were more than 5 sigma from the 
mean. A weighted mean and inverse variance were then calculated for all accepted pixels. 

The quality of the final images is displayed in Figures 4a-4d where we compare cutouts of our 
data in the HDF-N proper region with the Hubble Space Telescope images of Williams et al. (1996) 
The area with HK' data is shown in Figure 5. 

2.4. Object Detection and Measurement 

The catalog of objects was detected in the R and z' band images using Sextractor (Bertin & 
Arnouts 1996). Objects with three or more pixels rising more than 1.5a above the background 
were analyzed, but only those with 5a measurements of their aperture flux were output to the 
catalog. The significance of a detection was calculated using the RMS map generated during the 
image combination process. The absolute scaling of the RMS map was calculated by laying down 
random blank apertures away from objects on the images. This was done since the pixel-pixel noise 
underestimates the image noise. This method may somewhat overestimate the true noise since it 
also includes variations due to faint objects in the blank apertures. The area around bright stars 
was masked out to prevent erroneous detections and measurements from scattered light. 

We calculated both aperture and isophotal magnitudes for our images. Aperture magnitudes 
were selected as the primary magnitude to avoid biases introduced by galaxy morphology in isopho- 
tal magnitudes. To avoid filter dependent aperture corrections the images were smoothed with a 
Gaussian to match the median ratio of the total to the aperture magnitude in the worst seeing 
{U band) image. We chose this method because it ensures the apertures are sampling the same 
percentage of the total light in all bands and does not assume the PSF is gaussian. PSF matching 
provided a lower signal to noise because it includes the non-gaussian components of all seven PSF's 
which degraded the image quality. The aperture magnitudes were calculated on these smoothed 
images while the isophotal magnitudes were calculated on the un-smoothed images. A 3" diameter 
aperture was used for the photometry. An annulus with an inner radius of 6" and outer radius 
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of 12" was used for sky subtraction. The 5 sigma limiting magnitude cut was calculated using 
aperture magnitudes on the un-smoothed image to avoid losing faint objects. 

For the V band data we included archival data to increase the depth (Iwata et al. 2003). 
These data had better seeing (0.71" compared to 1.18" ); however they did not cover our entire 
field. Furthermore these data were taken with extremely long exposures which become non-linear 
at ~22nd magnitude To avoid position dependent image quality and saturation, we reduced the 
archival data separately and scaled them to have matching photometry. Where archival data existed 
and was unsaturated, it was used for the shape measurements. Photometry was performed on the 
two images separately and the fluxes combined using a weighted mean. 

2.5. Absolute Photometry 

The HDF-N is a heavily observed area of the sky with extremely accurate photometry in many 
bands. Many objects have known redshifts and spectral energy distributions making it easy to 
check the accuracy of our photometry. Furthermore it is difficult to accurately measure standard 
stars with Suprime-Cam because of the size of the telescope and overheads involved. Obtaining 
photometry on other telescopes would introduce problems similar to those in using the existing 
photometry. As such we have chosen not to use standard stars, but rather calibrate the data using 
the existing photometry. 

To obtain accurate absolute photometry in the AB system we matched our photometry to 
the HDF-N photometry of Fernandez-Soto et al. (1999) who carefully corrected for many subtle 
photometric effects in the HDF-N data. We converted the fluxes in the Fernandez-Soto et al. catalog 
to our filter system by linearly interpolating the fiux between the HDF-N filters. The effective 
wavelength of the filters were used for the interpolation. This method weights by transmission so 
any overlap between the Fernandez-Soto et al. (1999) filters and ours should be accounted for. We 
then compared our isophotal fluxes to the ones we interpolated and scaled them to match. Since the 
Fernandez-Soto et. al. fluxes were measured in an isophotal aperture, they vary systematically in 
magnitude with respect to the flux measured in a fixed aperture. To prevent this from introducing 
filter dependent offsets in our aperture fluxes we used the same magnitude range of 22-24 AB^ag in 
the U , B, V, R, I, and z' bands when comparing our aperture photometry to the Fernandez-Soto 
et. al. photometry. In the HK' band we were forced to use a range of 22-23.5 because of the 
shallower data. This may have introduced a small bias in our photometry for this band. We then 
scaled the aperture magnitudes to total magnitudes by comparing our isophotal and aperture fluxes 
for point sources in the / band. / band was chosen because it had the best image quality of the 
red bands and closely approximated the F814VF filter used in the HDF-N. To quantify the error 
in our zero point determination we measured the RMS scatter between our isophotal magnitudes 
and those of Fernandez-Soto et. al. The RMS scatter is 0.19 magnitudes which corresponds to and 
average error of 0.03 magnitudes in our photometric zero points exact numbers for all bands are 
given in table 3. Figures 6 and 7 compares our final photometry to the interpolated photometry of 
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Fernandez-Soto et al. (1999). 

Since an extensive catalog of spectroscopic redshifts (Cohen et al. 2000) exists for this field 
it is possible to compare our photometry to galaxy templates. This comparison was done using 
the photometric redshift code of Bcnftez (2000) which includes a calibration mode. First a best fit 
template at the known redshift was found for each object. The templates of Coleman et al. (1980) 
and Kinney et al. (1996) were used. Twenty intermediate spectra were interpolated between the 
main spectral templates to improve the fits. Then a weighed mean offset was calculated for each 
filter by comparing our photometry to the template photometry. These offsets were then applied to 
our photometry and the process repeated until the offsets converged. We did not apply offsets the B 
and / magnitudes because they are close to the HDF-N F450W and F814W filters respectively and 
hence we expect them to be correctly calibrated. The resulting offsets are small and are probably 
due to inaccuracies in the templates or filter profiles (see Table 3) . The largest offsets occur in the 
R, z', and HK' bands which are also those least like the HDF-N filters. The large offset in HK' 
may also be due to the different magnitude range used when scaling the aperture fluxes. 



2.6. Catalog Description 

Due to the large amount of data we have broken the catalog into several files which are available 
on the Internet-^^. There are two catalogs, one selected in the R band and the other in the z' band. 
Objects with an aperture magnitude error smaller than 5 sigma on the un-smoothed images in the 
selecting band were included. If an object was present in both catalogs it was removed from the z' 
catalog. Each object has a unique ID within each catalog, but both the R and z' catalog begin at 
1. Our data covers up to 0.4 square degrees in some bands, however around the edges the coverage 
is un-even, the astrometry may be distorted (see §2.2), and the detections become unreliable due 
to cosmic rays, reflections and other defects which could not be removed. As such we avoided these 
regions for our scientific objectives. All results in this paper are from the central 0.2 square degrees 
where the data is of equal depth in the U, B, V, R, I, and z' bands. Since others may find it 
useful we have included data for the entire field in supplementary files with the same formatting 
and content as the main catalog. However we strongly recommend using only the main R and z' 
catalog containing objects in the central 0.2 square degrees. No attempt was made to ensure the 
integrity of the data outside the central area. 

Each catalog consists of four files, the contents of which are listed in tables 4,5,6,7. The 
shape file contains information on the position and morphology of the object output by Scxtractor. 
The flag file contains flags for saturation, overlapping objects, and questionable or bad detections. 
The bad flag is set for objects with a Full Width at Half Max (FWHM) of or with abnormal 
magnitudes in the selection band. This flag was intended to mark questionable detections, however 
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some faint point sources may also have been flagged. The magnitude file contains aperture and 
isophotal magnitudes along with errors in the AB system for all bands. We used the Sextractor 
convention of -99 for an object which was outside the field or which could not be measured for 
other reasons. Badly saturated objects were also given a magnitude of -99. Objects with negative 
fluxes were assigned negative magnitudes with the absolute value of the magnitudes corresponding 
to the absolute value of the flux. The flux catalog contains aperture and isophotal fluxes along 
with errors and background levels. All measurements in the flux flle are in nano-Janskys (nJy). 

The area around bright objects has been removed due to false detections and problematic 
photometry in these areas. This was done by cutting out circles around USNO-A2.0 objects with 
magnitudes brighter than 15.5. A list of these regions and the size of each cutout is provided along 
with the catalog. 

3. Results 

3.1. Number Counts 

To determine number counts we constructed separate catalogs detected in each color using 
the un-smoothed images. Wc measured fluxes in 3" diameter apertures and applied an aperture 
correction to each band. We did not use the Sextractor best or isophotal magnitudes because they 
often act in a non-linear way at low signal to noise. The correction was calculated by measuring 
the median offset between Sextractor best magnitudes (Bcrtin &: Arnouts 1996), which estimate 
total magnitudes, and the aperture magnitudes for objects brighter than 24th magnitude (see Table 
3). Using aperture magnitudes may introduce a magnitude dependent bias in the number counts. 
However we saw no evidence that the median correction was changing with magnitude. Due to the 
difficulty in identifying stars at faint magnitudes we have provided both raw counts and counts with 
stars removed. Objects were identified as stars if they had a value of 0.9 or greater in the Sextractor 
star galaxy-separator. We did not remove stars in the HK' band because the star galaxy separator 
was unreliable in that band due to variable seeing across the image. The number counts shown in 
Figures 8-14 were normalized to a Euclidean slope as described in Yasuda et al. (2001). They agree 
with those reported by other authors (Metcalfe et al. 2001; McCracken et al. 2001; Huang et al. 
2001; Yasuda et al. 2001; Saracco et al. 2001; Djorgovski et al. 1995; Soifcr et al. 1994; Gardner 
et al. 1993). The Sloan Digital Sky Survey counts (Yasuda et al. 2001) were measured in n', g\ r', 
i', and z' which are close to our bands with the exception of g' and r' . The Sloan B band counts 
were generated by extrapolating between the u' , g\ and r' bands (Yasuda et al. 2001). To facilitate 
comparisons by future authors we have flt an exponential of the form N = 5iO"^("*"ff™*"''^) to our 
bright end counts where N is in number per square degree per magnitude. The results of these flts 
are quoted in table 8. The number count data are also provided in tables 9 and 10. 
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3.2. Selection of High Redshift Objects 

High redshift galaxies have unique colors due to the Lyman break at 912A and Lyman alpha 
absorption blueward of 1216A. We can select theses galaxies by choosing three filters which fall 
blueward of the Lyman break, between the Lyman break and Lyman alpha line, and redward of 
Lyman alpha respectively. Steidel ct al. (1999, 1995) successfully uses a set of customized filters to 
select galaxies at z ~ 3 and z ~ 4. Their spectroscopic foUowup is over 80% success in identifying 
z ~ 3 galaxies. Metcalfe et al. (2001) has also used these criteria on the Herschel Deep Field and 
HDF-N. 

The number and type of object selected are strongly dependent on the filters and selection 
criteria used. Changing these parameters can severely effect the success in identifying high redshift 
galaxies and change the redshift range selected. For comparison to the previous work we need to 
select the same population as Steidel et al. (1999, 1995) and Metcalfe et al. (2001). However our 
band passes differ significantly from theirs. We corrected this by calculating the expected colors 
for the galaxies they have selected. These colors are calculated by integrating the galaxy templates 
of Coleman et al. (1980) and Kinney ct al. (1996) moderated by our filter response profile. We 
correct these templates for Ly^ and Ly^ absorption following the prescription of Madau (1995). 
We also calculate the expected colors of stars using the spectral library of Pickles (1998). We then 
set a selection criteria which will select similar galaxies to the previous work but avoids stars. The 
resulting selection is shown in Figures 17 and 19 along with color-color tracks for various galaxy 
types. The selection is also shown along with our data in figures 16 and 18. 

The surface density we measure using this selection is similar to that measured by other authors 
(Metcalfe et al. 2001; Steidel et al. 1999) (see Figures 22,23 and tables 11,12. Due to our redder 
U band filter the z ~ 3 selection suffers more contamination from stars and low redshift galaxies 
than that of Steidel et al. (1999, 1995). We attempted to remove stars by selecting objects with 
R < 23.5 and Rfwhm > 1-3" as describe in McCracken et al. (2001), however the contamination 
from galaxies remains. At magnitudes fainter than i? = 23 our measurements are within the one 
sigma error bars reported by Steidel et al. (1999, 1995) and Metcalfe et al. (2001). Our B band 
counts are consistent with but systematically lower than those of Steidel et al. (1999) or Metcalfe 
et al. (2001) which is likely due to cosmic variance. 

The selection at redshifts higher than 2; ~ 4 is significantly more difficult. The contamination 
for z ~ 3 galaxies comes from low redshift galaxies where the co-moving volume is small. For 
z > 5 galaxies the contamination comes from z ~ 1 where the effective co-moving volume is very 
large. Several authors have recently reported the discovery of large numbers of z > 5 galaxies using 
color-color selection (Iwata et al. 2003; Lehnert & Bremer 2003). However we find their two color 
selection yields a large number of lower redshift galaxies. Figure 15 shows an example for a redshift 
5.5 galaxy. Elliptical galaxies at 1 < 2; < 1.5 have similar colors to a z > 5.5 galaxy in the R, 
I, and z' bands. The only clear difference is in the bluer or redder bands. In particular the use 
of an (R — I) selection will include many low redshift galaxies. A similar problem will occur at 
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redshift 6 using an (/ — Z) selection. There is no easy solution, either one must go more than three 
magnitudes deeper in the bluest band than the reddest, or obtain deep near IR imaging, both of 
which are time intensive for large surveys. 

Songaila et al. (1990) sets out a clear criteria for selecting high redshift galaxies based on the 
Lyman alpha absorption shown in equation 2. 



Using these criteria and our multi-color data we can test the findings of other authors. To test 
the level of contamination we have restricted ourselves to z' > 25 which is bright enough for the 
expected contaminating sources to be detected at 2a in U or B bands. 

Iwata et al. (2003) imaged the same area with Suprime-Cam in the same V, I, and z' filter 
we used. However they calibrated their data in a different manner. We adopted {V — I) > 2.4 and 
{V — I) > 7{I — z') — 0.2 as our selection criteria. This avoids selecting z ~ 1 galaxies and late 
type stars. We find that our counts arc consistent with, but lower than the values of Iwata ct al. 
(2003) (see Figure 24 and Tables 13 and 14). Moving our selection slightly redder in (/ — Z) color 
increases the number of objects selected, making our numbers more consistent with Iwata et al. 
(2003). However all the additional objects are detected inU oi B band. Even when objects detected 
in U OT B are removed, many of the remaining objects in the ACS GOODS field morphologically 
appear to be stars (this will be quantified elsewhere). These results suggest that a combination of 
z ~ 1 galaxies and late type stars are contaminating the z > 5.5 selection. 

Figure 21 shows that the selection criteria of Lehnert & Bremer (2003) is very likely selecting 
low redshift galaxies or stars. Wc find that 95% of the objects selecting using this criteria are 
detected at 2a in V. This selection is particularly problematic since it may include star forming 
galaxies at z < 2 where [O III] or [O II] can be mistaken for Lyman a. Finally the (I — Z) > 1.5 
selection used by Stanway et al. (2003), selects objects which are detected at 2a in V band 53% 
of the time. In contrast only 14% of objects selecting using the {B — R), {R — I) selection were 
detected in U band. Based on our number counts the probability of an object at 28*'* magnitude 
falling into a 3" aperture by chance is 37% in V band and 4% in U band. This implies that color 
selection may be overestimating the number of high redshift galaxies and hence the star formation 
rates at z > 5. Wc shall give a more extensive discussion of this issue with photometric redshift 
estimates in a future paper (Capak et al. 2003). 




(2) 



For z > 5 galaxies it requires (V^ — /) > 2.4 and no detection blueward of the 912A break. 
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3.3. Conclusion 

We have compiled a deep, multi-color catalog over 0.2 square degrees in the HDF-N region. 
These data will provide an invaluable basis for understanding the formation and evolution of galaxies 
by providing a large sample of galaxies which can be studied in various ways. The raw number 

counts and counts of z ~ 3, and z ~ 4 galaxies agree with those of other authors. However the 
selection of z > 5 galaxies suffers from significant contamination by low redshift galaxies. The 
current deep blue imaging is required to reliably selected these galaxies and hence estimate the star 
formation rate. 

We'd like to thank Nick Kaiser for providing his software and supporting it. We'd also like 
to thank H. McCracken and N. Metcalfe for providing their number count data. The work of C. 
McNally and D. Stern was carried out at Jet Propulsion Laboratory, California Institute of Tech- 
nology, under a contract with NASA. This work was supported by NASA grant DF1-2001X(L.L.C) 
and NSF grants ASTOO-71208 (E.M.H), AST 00-84816 (L.L.C.), and AST 99-83783 (A.J.B.). 
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Fig. 1. — Filter profiles multiplied by the quantum efficiency of the detectors used from our survey 
(top) and the HDF-N proper (Fernandez- Soto et al. 1999) (bottom). All profiles were normalized 
to have a peak transmission of one. 
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Scatter with respect to radio positions 



Astrometric residuals for radio sources 
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Fig. 2. — Astrometric residuals with respect to the radio catalog of Richards (2000). The inner box 
on the left hand side marks the central 0.2 square degrees in the residual map. The vectors point 
in the direction of the residual and have and amplitude of 60 times the residual in arc seconds. 
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Fig. 3. — Astrometric residuals with respect to the USNO-Bl.O (Monet et al. 2003). The inner box 
on the left hand side marks the central 0.2 square degrees in the residual map. The vectors point 
in the direction of the residual and have and amplitude of 60 times the residual in arc seconds. 
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Fig. 4a. — On the left is a cutout image of our data around the HDF-N proper in U band, and on 
the right is the Hubble space telescope band image for comparison (Williams et al. 1996). 

Note the depth and quality of the ground based image. 
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Fig. 4b. — On the left is a cutout image of our data around the HDF-N proper in B band, and on 
the right is the Hubble space telescope F450W band image for comparison (Wilhams et al. 1996). 
Note the depth and quahty of the ground based image. 
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Fig. 4c. — On the left is a cutout image of our data around the HDF-N proper in V band, and on 
the right is the Hubble space telescope F606W band image for comparison (Williams et al. 1996). 
Note the depth and quality of the ground based image. 
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Fig. 4d. — On the left is a cutout image of our data around the HDF-N proper in / band, and on 
the right is the Hubble space telescope FSIAW band image for comparison (Williams et al. 1996). 
Note the depth and quality of the ground based image. 
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Hawaii-HDF-N AB magnitude 

Fig. 6. — Comparison of our photometry to that of Fernandez-Soto et al. (1999) from the HDF-N. 
The magnitudes have been converted to our filter system by linearly interpolating the flux between 
the HDF-N filters. The RMS scatter between our isophotal magnitudes and those of Fernandez- 
Soto et al. is 0.19 magnitudes on average which corresponds to an average error of 0.03 magnitudes 
in the zero points. 
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Fig. 7. — Comparison of our photometry to that of Fernandez-Soto et al. (1999) from the HDF-N. 
The magnitudes have been converted to our filter system by linearly interpolating the flux between 
the HDF-N fllters. The RMS scatter between our isophotal magnitudes and those of Fernandez- 
Soto et al. is 0.19 magnitudes on average which corresponds to an average error of 0.03 magnitudes 
in the zero points. 
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Fig. 8. — U band number counts from this work, the Herschel Deep Field (Metcalfe et al. 2001), 
the HDF-N (Metcalfe et al. 2001), and the Sloan Digital Sky Survey (SDSS) (Yasuda et al. 2001). 
Objects were identified as stars if they had a value of 0.9 or greater in the Sextractor star galaxy- 
separator. 
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Fig. 9. — B band number counts from this work, the Calar Alto Deep Imaging Survey (Huang et al. 
2001), the Canada Prance Deep Survey (CFDS)(McCracken et al. 2001), the HDF-N (Metcalfe et al. 
2001), and the Sloan Digital Sky Survey (SDSS) (Yasuda et al. 2001). Objects were identified as 
stars if they had a value of 0.9 or greater in the Sextractor star galaxy-separator. 
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Fig. 10. — V band number counts from this work, the HDF-N (Metcalfe et al. 2001), and the Sloan 
Digital Sky Survey (SDSS) (Yasuda et al. 2001). 
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Fig. 11. — R band number counts from this work, the Calar Alto Deep Imaging Survey (Huang 
et al. 2001), the Herschel Deep Field (Metcalfe et al. 2001), the HDF-N (Metcalfe et al. 2001), and 
the Sloan Digital Sky Survey (SDSS) (Yasuda et al. 2001). Objects were identified as stars if they 
had a value of 0.9 or greater in the Sextractor star galaxy-separator. 
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Fig. 12. — I band number counts from this work, the Canada France Deep Survey 
(CFDS)(McCracken et al. 2001), the Herschel Deep Field (Metcalfe et al. 2001), the HDF-N (Met- 
calfe et al. 2001), and the Sloan Digital Sky Survey (SDSS) (Yasuda et al. 2001). Objects were 
identified as stars if they had a value of 0.9 or greater in the Sextractor star galaxy-separator. 
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Fig. 13. — z' band number counts from this work and the Sloan Digital Sky Survey (SDSS) (Yasuda 
et al. 2001). Objects were identified as stars if they had a value of 0.9 or greater in the Sextractor 
star galaxy-separator. 
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Fig. 14. — K band number counts from this work, the Chandra Deep Field South (CDF-S) (Saracco 
et al. 2001), the HDF-N (Saracco et al. 2001), and several other authors (Bershady et al. 1998; 
Moustakas et al. 1997; Djorgovski et al. 1995; Soifer et al. 1994; Gardner et al. 1993) are over 
plotted. An offset of 1.859 was applied to Vega magnitudes to convert them to the AB system. 
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Fig. 15. — The expected Spectral Energy Distribution (SED) of a z = 5.5 star forming, z = 1.5 
Elliptical, and z = 1.0 Elliptical galaxy from the Coleman et al. (1980) library are shown along 
with an M9 star from Pickles (1998). These SED's were corrected for intergalactic absorption using 
Madau (1995). The filter profiles of our survey are overlaid for comparison. Note the similarity 
of the SED's in R, I, and z' bands. Note that differentiating the z = 5.5 Irr from the z = 1.0 or 
z = 1.5 elliptical requires multi-color information and is most sensitive in the C/, -B, and near-IR 
bands. The maximum information comes from the largest wavelength range sample. 
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Fig. 16. {D — R) vs. {U — B) color-color plot for all objects with Bab > 26. The region 
indicated by heavy lines in the top right is where we expect to find z ~ 3 objects similar to those 
of Steidel et al. (1995). The light solid lines show the expected evolution of Coleman et al. (1980) 
galaxy templates and the dotted lines the Kinney et al. (1996) SB2 and SB3 galaxy templates. 
The expected colors of stars from Pickles (1998) are over plotted as stars. The downward facing 
triangles at {U — B) = 6 are objects which were not detected in U. 
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Fig. 17. — (-B — R) vs. {U — B) color-color plot with galaxy tracks. The region indicated by heavy 
hnes in the top right is where we expect to find z ~ 3 objects similar to those of Steidel et al. (1995). 
The light solid lines show the expected evolution of Coleman et al. (1980) galaxy templates and 
the dotted lines show the Kinney et al. (1996) SB2 and SB3 galaxy templates. The expected colors 
of stars from Pickles (1998) are over plotted as stars. Circles indicate galaxies in the 2.5 < z < 3.5 
redshift range selected by Steidel et al. (1995). 



-37- 




Fig. 18. — {R — I) vs. {B — R) color-color plot for all objects with Rab > 26. The region indicated 
by heavy lines in the top left is where we expect to find z ~ 4 objects similar to those of Steidel 
et al. (1999). The light solid lines show the expected evolution of Coleman et al. (1980) galaxy 
templates and the dotted lines show the Kinney et al. (1996) SB2 and SB3 galaxy templates. The 
expected colors of stars from Pickles (1998) are over plotted as stars. The downward facing triangles 
at {B — R) = 6 are objects which were not detected in B. 
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Fig. 19. — (R — I) vs. (B — R) color-color plot with galaxy tracks. The region indicated by heavy 
hnes in the top left is where we expect to find z ~ 4 objects similar to tliose of Steidel et al. (1999). 
The light solid lines show the expected evolution of Coleman et al. (1980) galaxy templates and 
the dotted lines show the Kinney et al. (1996) SB2 and SB3 galaxy templates. The expected colors 
of stars from Pickles (1998) are over plotted as stars. Circles indicate galaxies in the 3.5 < z < 4.8 
redshift range selected by Steidel et al. (1995). 
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Fig. 20. — (/ — z') vs. {V — I) color-color plot for all detected objects. The region indicated by 

heavy lines in the top left is where we expect to find z ~ 5.5 objects based on the criteria of Iwata 
et al. (2003). Wc have moved the cuts to {V - I) >= 2.4 and {V - I) > 7(1 - Z) - 0.2 to avoid 
contamination from lower redshift galaxies and late type stars as well as account for differences in 
the photometric system. The expected colors of stars from Pickles (1998) are over plotted as stars. 
The downward facing triangles dX {V — I) = 6 are objects which were not detected in V. 
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Fig. 21. — {I—z') vs. (R — I) color-color plot for all detected objects. The heavy line at (/? — /) = 1.5 
indicates the selection criteria of Lehnert & Bremer (2003) and the dashed line at (I — Z) = 1.5 is 
the selection criteria of Stanway et al. (2003). The light solid lines show the expected evolution of 
Coleman et al. (1980) galaxy templates and the Kinney et al. (1996) SB2 and SB3 galaxy templates. 
The expected colors of stars from Pickles (1998) are over plotted as stars. The downward pointing 
triangles at (i? — /) = 6 are objects which were not detected in R and leftward pointing triangles 
at (/ — Z) = 6 are objects which were not detected in I. 
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Fig. 22. — The raw sky density of U band dropouts uncorrected for contamination or incompleteness 
are shown. Error bars are la poison fluctuations for Metcalfe et al. (2001) and our data. The errors 
for Steidel et al. (1999) include an estimate of cosmic variance from their multiple fields. Much of 
the variation between the counts is likely due to cosmic variance. 
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Number density of B dropouts 
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Fig. 23. — The raw sky density of B band dropouts uncorrected for contamination or incompleteness 
are shown. Error bars are la poison fluctuations for Metcalfe et al. (2001) and our data. The errors 
for Steidel et al. (1999) include an estimate of cosmic variance from their multiple fields. Much of 
the variation between the counts is likely due to cosmic variance. 
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Number density of V dropouts 
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Fig. 24. — The raw sky density of V band dropouts and an estimate of the density corrected for 
contamination are shown for the Hawaii-HDF-N. Neither measurement is corrected for incomplete- 
ness. Error bars are la poison fluctuations. The counts of Iwata et al. (2003) are higher due to a 
higher rate of contamination. 



Table 1. Telescopes used, Exposure Times, and Dates of Observations 



Band 


Integration^ 


Telescope 


Typical 


Date oi 




iime (h) 


Used 


Integration (s) 


Observations 


u 


28.5 


KPNO 4m 


1800 


2002, March, 9-13 


B 


1.7 


Subaru 8.3m 


600 


2001, February, 27,28 


V 


6.4 


Subaru 8.3m 


1200 


2001, February, 23,24 




2.0 




600 


2001, March 21-23 


R 


5.2 


Subaru 8.3m 


480 


2001, February, 27,28 








600 


2001, March, 21-23 


I 


2.9 


Subaru 8.3m 


300 


2001, February, 23,24 








300 


2001, March, 21-23 








300 


2002, April, 5,11-14 


z< 


3.9 


Subaru 8.3m 


180 


2001, February, 23,24,27,28 








240 


2000, March, 21-23 


HK' 


0.43 


UH 2.2m 


120 


1999-2002 



^This is the average integration time for a pixel. This varies significantly in the HK' image. 



Table 2. Data Quality, Depth, and Coverage 



Band 


Seeing 


5a limit'^ 


Total area 


Deep area 




Date of 




arcsecond 


(AB mag) 


Covered (deg^) 


Covered (deg^) 




Observations 


T T 

U 


1.2b 


27.1 


0.40 


u.ob 


2002, 


March, 9-13 


B 


U.71 


26.9 


0.27 


0.20 


2001, 


l^ebruary, 27,28 


V 


0.71 


26.8 


0.27 


0.20 


2001, 


February, 23,24 




1.18 




0.40 


0.20 


2001, 


March 21-23 


R 


1.11 


26.6 


0.27 


0.20 


2001, 


February, 27,28 








0.40 


0.20 


2001, 


March, 21-23 


I 


0.72 


25.6 


0.27 


0.20 


2001, 


February, 23,24 








0.40 


0.20 


2001, 


March, 21-23 








0.40 


0.20 


2002, 


April, 5,11-14 


z' 


0.67 


25.4 


0.27 


0.20 


2001, 


February, 23,24,27,28 








0.40 


0.20 


2000, 


March, 21-23 


HK' 


0.87 


22.1 ^ 


0.11 


0.11 


1999-2002 



^Measured by randomly placing 10'^ 3" diameter apertures on an image with detected objects 
masked out. The apertures were placed at least 3" away from the positions of detected object. 

^This is the mode of the depth across the field. In a 9'x 9' area around the HDF-N the 5a limit 
is 22.8. 



Table 3. Photometric Information 



Band 


Zeropoint 


Offsets from 


Aperture 


Offsets to 


Central 


Bandpass A 




Error^ 


SED fitting; 


(norrpption ^ 


Vpfa svstpm^ 


WavelpTiP'tli A 




U 


0.031 


-0.008 


-0.311 


0.719 


3647.65 


387.16 


B 


0.028 


0.000 


-0.159 


-0.077 


4427.60 


622.05 


V 


0.028 


-0.005 


-0.282 


0.023 


5471.22 


577.36 


R 


0.026 


0.061 


-0.262 


0.228 


6534.16 


676.39 


I 


0.016 


0.000 


-0.126 


0.453 


7975.89 


792.89 


z' 


0.018 


-0.065 


-0.147 


0.532 


9069.21 


802.63 


HK' 


0.042 


-0.181 


-0.233 


1.595 


18947.38 


5406.22 



'^Includes the estimated error of 0.01 magnitudes in the Fernandez-Soto et al. (1999) catalog. 
^Used in the calculation of number counts. 

^Calculated by integrating an AO star SED multiplied by our filter profile. 



Table 4. Contents of Shape File 



Column Number Shape Catalog Item 



1 


ID 


2 


RA(J2000) 


3 


DEC(J2000) 


4 


X position on image 


5 


y position on image 


6 


U FWHM 


7 


U Semi-Major Axis ^ 


8 


U Semi-Minor Axis 


9 


B FWHM 


10 


B Semi-Major Axis ^ 


11 


B Semi-Minor Axis ^ 


12 


V FWHM <^ 


13 


V Semi-Major Axis ^ 


14 


V Semi-Minor Axis 


15 


R FWHM ^ 


16 


R Semi-Major Axis ^ 


17 


R Semi-Minor Axis * 


18 


I FWHM ^ 


19 


I Semi-Major Axis ^ 


20 


/ Semi-Minor Axis ^ 


21 


z' FWHM * 


22 


z' Semi-Major Axis * 


23 


z' Semi-Minor Axis ^ 


24 


HK' FWHM ^ 


25 


HK' Sciiii-Major Axis'"" 


20 


/././\' S(uiii-.Miii(ir .Axis'^ 



^This is the output from the Scxtractor FWHM-IMAGE measurement. The FWHM is calculated in a different manor than 
the Semi-Major and Semi-Minor axis which which are outputs from the A-IMAGE and B-IMAGE measurements. All three of 
these measurements were made on the un-smoothed images with respect to centroids in the detection images and have been 
scaled to arcseconds. A value of indicates that the measurement could not be determined. For more information on these 
parameters see Bertin & Arnouts (1996). 
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Tablc 5. Contents of Flag File 



( '(.iliuiiii Xuiiil)('i' 


Flag Calalog Item 


( 'oillUKUil S 


1 


ID 




2 


Bad 


Set to 1 if the FWHM is in the detecting band 


3 


U Saturated 


Set to 1 if any pixel is above saturation the limit 


4 


B Saturated* 


Set to 1 if B < 20.49 


5 


V Saturated'* 


Set to 1 if y < 19.85 


6 


_R Saturated^ 


Set to 1 if K < 20.25 


7 


I Saturated* 


Set to 1 if / < 20.15 


8 


z' Saturated* 


Set to 1 if Z < 19.45 


9 


HK' Saturated 


Set to 1 if any pixel is above saturation the limit 


10 


N Overlapping 


Number of objects detected within 3" of this object 



^Saturated pixels behave in a very non-linear fashion on Suprime-Cam and can drop in value after reaching the saturation 
limit. We adopted these magnitude cuts by measuring the magnitude at which point sources became non-Gaussian. 
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Table 6. Contents of Magnitude File 



Column Number 


Magnitude Catalog Item'' 


1 


ID 


2 


U Aperture Magnitude 


3 


T T \ J_ TIT "j-l 1 — ^ 

U Aperture Magnitude trror 


4 


(7 Isopnotal Magnitude 


5 


U Isophotal Magnitude Error 


6 


B Aperture Magnitude 


7 


B Aperture Magnitude Error 


8 


B Isophotal Magnitude 


9 


B Isophotal Magnitude Error 


10 


V Aperture Magnitude 


11 


V Aperture Magnitude Error 


12 


V Isophotal Magnitude 


13 


V Isophotal Magnitude Error 


14 


R Aperture Magnitude 


15 


R Aperture Magnitude Error 


16 


R Isophotal Magnitude 


17 


R Isophotal Magnitude Error 


18 


I Aperture Magnitude 


19 


/ Aperture Magnitude Error 


20 


/ Isophotal Magnitude 


21 


/ Isophotal Magnitude Error 


22 


z' Aperture Magnitude 


23 


z' Aperture Magnitude Error 


24 


z' Isophotal Magnitude 


25 


z' Isophotal Magnitude Error 


26 


HK' Aperture Magnitude 


27 


HK' Aperture Magnitude Error 


28 


HK' Isophotal Magnitude 


29 


HK' Isophotal Magnitude Error 
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Table 7. Contents of Flux File 



Column Number Flux Catalog Item 



1 


ID 


2 




3 




4 


// A I'lOT'tii fo r'k'CTV*"!] 1 nrl 


5 


t/ Isopliota,! Flux 


g 


TJ T^nnVintal Flii'v Terror 


7 


f/ Isophot&l Backgroimd 


g 




g 


R A T^f y+llT'l^ TillTV niTTlTr 


10 


S Aperture Bctckj^round 


11 


R TtionVinfa] T^lii">f 


12 


R Tsjnnlnntal Flmr PiTTrtr 


13 




14 


T/^ Anprtnrp "Flnv 


15 


1/ A "nr>T'("ii fo H^liiY rOrmv 

V JA. LJt_.i L Lii X 1 \Jlj\. XL;i i LJX 


16 


1/ A"^^/^T^"11T^^ \-\ OicA^ (vyrw i ti *H 

V Jt. Ul^i Lr Lii t, J_)CiL,iS.^l tJ LlliU. 


17 


Isophotal Plux 


18 


Isopliotcil Flux Error 


19 




20 


X L r^vj^i. \j i_ix ^ X X LXA. 


21 




22 


R Aperture Bcickground 


23 


Isopliotcil Flux 


24 


R Term n rtt" n 1 Ti liiv Hiwi"!!" 


25 


R Isophotal Background 


26 


/ Aperture Flux 


27 


/ Aperture Flux Error 


28 


I Aperture Background 


29 


/ Isophotal Flux 


30 


/ Isophotal Flux Error 


31 


I Isophotal Background 


32 


z' Aperture Flux 


33 


z' Aperture Flux Error 


34 


z' Aperture Background 


35 


z' Isophotal Flux 


36 


z' Isophotal Flux Error 


37 


z' Isophotal Background 


38 


HK' Aperture Flux 


39 


HK' Aperture Flux Error 


40 


HK' Aperture Background 


41 


HK' Isophotal Flux 


42 


HK' Isophotal Flux Error 




/././v' ls(>i)li()l aJ Ljackgrouiid 



''AH flux values are in nano-Janskys (njy). 
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Table 8. Fit of the form N = filO^^^^^ magnitude) ^j^g number counts^ 



Band 


Fit Range in AB magnitudes 


A 


Error in A 


logio{B) 


Error in logio{B) 


U 


20.0-24.5 


0.526 


0.017 


-8.61 


0.37 


B 


20.0-25.5 


0.450 


0.008 


-6.67 


0.18 


V 


20.0-25.0 


0.402 


0.004 


-5.50 


0.09 


R 


20.0-25.0 


0.361 


0.004 


-4.36 


0.08 


I 


20.0-25.0 


0.331 


0.008 


-3.54 


0.17 


z' 


20.0-25.0 


0.309 


0.006 


-2.98 


0.13 


HK' 


16.0-19.0 


0.378 


0.037 


-4.15 


0.65 


HK' 


19.0-22.0 


0.393 


0.019 


-4.31 


0.39 



has units of number degree ^ 0.5 mag 



Table 9. Raw number counts in N degree ^ 0.5 mag 



Magnitude 


U 


B 


V 


R 


/ 


z' 


HK' 




41 it 14 














18.75 


46 it 15 










454 + 48 


QftQ _L on 

ouf? ^ oy 




77 -t- 9n 








CI 1 -LCI 

ui-j- m oi 


ooo ZL oo 


1 f\7'\ -A- 1 9A 




QS -t- 99 




'iAf\ A- 49 




y±^ m oo 


1 m 9 _|_ 79 


zyyu jz ±ou 


20 25 


1 1T^ -t- 9*^ 
J.UO m zo 


oiD in 




QC9 _i_ ap. 


1 99Q _i_ 7Q 
j.zzy m 1 y 




/1«Q9 _|_ 9nR 


9n 7^ 

zu. / o 


9^*^ -t- '^fi 
zuo in ou 


^tjo in ou 


79*^ A- (\^ 
1 zo in oi- 


J- zoo m ou 


91 on 4- infi 
ziyu m Auu 


9snn -1- 1 9n 


7140 -1- 9c;7 


21.25 


371 ± 43 


712 it 60 


1100 it 75 


2190 it 106 


3259 it 129 


3797 ± 140 


12317 ± 337 


21.75 


557 it 53 


1224 it 79 


1777 it 95 


3048 ± 125 


4918 ± 159 


6044 ± 176 


15655 it 380 


22.25 


1048 ± 73 


2071 it 103 


2588 it 115 


4551 ± 153 


7113 ± 191 


8854 ± 213 




22.75 


2371 it 110 


3435 it 133 


4282 it 148 


6850 it 188 


10301 it 230 


12135 lb 250 




23.25 


4096 ± 145 


6468 it 182 


7067 it 191 


10585 it 233 


13984 it 268 


16624 it 293 




23.75 


8389 ± 208 


11887 it 247 


11484 it 243 


16268 it 289 


20380 it 324 


22622 lb 341 




24.25 


14284 it 271 


19647 it 318 


18727 it 311 


24224 it 353 


29488 it 390 


33802 lb 417 




24.75 


19982 ± 321 


30361 it 396 


28858 it 386 


35925 ± 430 


42827 it 470 


45235 lb 483 




25.25 


29824 it 392 


45064 ± 482 


40869 ± 459 


48474 ± 500 


51016 it 513 


44718 lb 480 




25.75 


42135 it 466 


60542 it 559 


44073 ± 477 


61410 it 563 








26.25 


51093 it 513 


70141 it 601 


27918 it 379 


51863 it 517 








26.75 


39361 it 450 


55831 it 537 












27.25 


21568 it 333 















Table 10. Number counts with stars removed in N degree ^ 0.5 mag 



Magnitude 


U 


B 


V 


R 


/ 


z' 




^ -t- ^ 












lO. / O 


o in o 










9*^9 -t- "iA 




O HI o 








'i'i^ -4- A^ 
ooo m ^1 


f^9i -t- ^^ 




o in o 








fi7S -t- "^4 
o / o HI 


1 CiA'i -t- 7'^ 


20.25 


5 ± 5 


154 it 28 


253 it 36 


625 it 56 


1203 it 78 


1839 ± 97 


20.75 


87 ± 21 


346 it 42 


464 it 49 


1198 it 78 


2076 it 103 


2898 it 122 


21.25 


191 ± 31 


557 it 53 


873 it 67 


2087 ± 103 


3104 ± 126 


4391 it 150 


21.75 


418 ± 46 


1079 it 74 


1415 it 85 


2794 ± 120 


4990 it 160 


6297 it 180 


22.25 


821 ± 65 


1932 ± 99 


2464 it 112 


4510 it 152 


7501 it 196 


9738 i 224 


22.75 


2014 ± 102 


3290 it 130 


3952 it 142 


7000 it 190 


10973 it 238 


12688 it 256 


23.25 


3699 ± 138 


6344 it 181 


6840 it 187 


10823 it 236 


15255 it 280 


18055 it 305 


23.75 


7728 ± 199 


11784 it 246 


11448 it 243 


16516 ± 292 


22617 ± 341 


24658 ± 356 


24.25 


13111 it 260 


19564 it 317 


18918 it 312 


24224 it 353 


32516 it 409 


36705 it 435 


24.75 


18737 ± 311 


30274 ± 395 


28977 ± 386 


34618 it 422 


46170 ± 488 


46485 it 490 


25.25 


28316 ± 382 


44956 it 481 


41376 it 462 


46289 it 489 


50778 it 512 




25.75 


41386 it 462 


60434 it 558 


43453 it 473 


59370 it 553 






26.25 


51062 it 513 


69986 it 601 
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Table 11. Number density of U dropouts 



R magnitude 


With Stars 




Without Stars 




N arcmin""^ 0.5 mag""*^ 


Error 


N arcmin ^ 0.5 mag ^ 


Error 




22.25 


0.017 


0.005 


0.003 


0.002 


22.75 


0.030 


0.006 


0.017 


0.005 


23.25 


0.040 


0.007 


0.031 


0.007 


23.75 


0.095 


0.012 


0.095 


0.012 


24.25 


0.214 


0.017 


0.214 


0.017 


24.75 


0.450 


0.025 


0.450 


0.025 


25.25 


0.661 


0.031 


0.661 


0.031 


25.75 


0.653 


0.030 


0.653 


0.030 
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Table 12. Number density of B dropouts 



/ magnitude 


N arcmin ^0.5 mag ^ 


Error 


23.25 


0.0 




23.75 


0.014 


0.004 


24.25 


0.035 


0.007 


24.75 


0.099 


0.011 


25.25 


0.198 


0.017 
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Table 13. Number density of V dropouts. 



/ maguitudo " 


N arcmin ^ 0.5 mag ^ 


Error 


23.75 


0.011 


0.004 


24.25 


0.019 


0.005 


24.75 


0.051 


0.008 


25.25 


0.059 


0.009 


25.75 


0.061 


0.009 



^Ic = I — 0.453 for comparison to Iwata et al. (2003) 



Table 14. Number density of V dropouts not detected at 2a in U or B with V — I > 2A 



I magnitude N arcmin ^ 0.5 mag ^ Error Percentage of Interlopers 



23.75 0.001 0.001 90.9 

24.25 0.010 0.004 47.4 

24.75 0.038 0.007 25.5 

25.25 0.047 0.008 20.3 



